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Abstract. We give a detailed analysis of the oscillation formula within the context of the 
wave packet formalism. Particular attention is made to insure flavor eigenstate creation 
in the physical cases (Ap 7^ 0). This requirement imposes non instantaneous particle 
creation in all frames. It is shown that the standard formula is not only exact when the 
mass wave packets have the same velocity, but it is a good approximation when minimal 
slippage occurs. For more general situations the oscillation formula contains additional 
arbitrary parameters, which allows for the unknown form of the wave packet envelope. 

PACS. 12.15.F - 14.60.Pq 

I. INTRODUCTION 

In the last few years, the great interest in neutrino physics, and in particular in neutrino masses, 
has stirred up an increasing number of theoretical works upon the quantum mechanics of oscillation 
phenomena [1-3]. Notwithstanding the exceptional ferment in this field, the situation is still confused, 
and the conceptual difficulties hidden in the oscillation formulas represent an intriguing, and some- 
times embarrassing, challenge for physicists. The most controversial point in discussing the quantum 
mechanics of particle oscillations is represented by the derivation of formulas containing extra factors 
in the oscillation length [4 12]. In this paper, we review the source of these factors and show why 
in the wave packet formalism with minimal slippage one finds, independently from the kinematical 
assumptions, the standard oscillation probability. However, in more general situations one does indeed 
find different oscillation formulas. This question is of course essential if one wants to derive consistent 
mass differences from oscillation phenomena. It may well be that different experiments, or even sets of 
measurements within a given experiment (e.g. atmospheric neutrino data), involve different oscillation 
formulas, and this must be taken into account. 

One of the basic assumptions in neutrino physics is that only flavor eigenstates are destroyed or 
created. Now, in the wave function formalism this is a problem which, in our opinion, has not yet 
been satisfactorily solved. The most common approach is to assume instantaneous creation within 
the context of an equal momentum hypothesis [13]. Unfortunately, there is no physical Lorentz frame 
in which this occurs. Some authors have even tried to bypass the problem by re-defining the flavor 
eigenstates according to convenience [14] . In this paper we shall describe how to achieve this obligatory 
condition within the wave packet formalism. It will automatically imply non instantaneous creation 
of the wave packet for any physical production process. 

In our analysis, we shall for simplicity work within a two flavor/mass mixing model. The structure 
of the paper is as follows: 
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• In Section II, we recall the arguments in favor of extra factors in the oscillation formula. We shall 
in particular note that an essential assumption in these derivations is that the flavor eigenstate is 
identical, including the phase, at each point of creation. 

• In Section III, we recall the equal momentum wave packet derivation of the standard oscillation 
formula. 

• We discuss the physical kinematics for particle creation in Section IV. 

• The extension from the equal momentum case to an arbitrary case is tackled, and solved, in Section 
V. In this Section, we also show that, while physically essential, this generalization does not in itself 
invalidate the standard oscillation formula. 

• In Section VI, we return to the question of non-standard oscillation formulas by showing the results 
of a multiple peak (specifically a two-peak) model with substantial slippage. 

• We draw our conclusions in Section VII. 



II. EXTRA FACTORS IN THE OSCILLATION FORMULAS 

In the quantum mechanics of particle oscillation, substantial mathematical simplification results from 
the assumption that the space dependence of the wave functions is one-dimensional, hence, in what 
follows, we shall use this simplification. Flavor oscillations are observed when a source creates a particle 
which is a mixture of two or more mass eigenstates. The main aspects of oscillation phenomena can 
be understood by studying the two flavor problem 

v a \ _ / cos 9 sin 9 \ ( v x 
Vp J \ -sin 9 cos 9 J \u 2 

where v a and u p are flavor eigenstates and v 1 and v 2 are mass eigenstates. 

Suppose we have a physical system whose initial state is represented by the flavor state v a . At 
later times, the probability to find the flavor state Vp is conventionally expressed in terms of the 
mixing angle 9 and of the relative phase A<f> by 

, ,„„1- cos A$ 
P (u a -> v„; t) = sm 2 29 . 

It is to be noted that the above formula ignores any possible effects involving the shape of the wave 
functions. The Lorentz invariant phase factor A<!> is usually given in terms of the distance L (travelled 
in the time T), of the mass difference Am 2 , and of the energies mean value E. Due to the relativistic 
nature of neutrinos, this phase is evaluated by considering T«L and p 1>2 ~ E 1-2 , i.e. 

A$ = T AE -LApwL (AE — Ap) ps L Am 2 /2E . (1) 

By using this phase difference, one gets the well-known expression [15-19] 

P {u a -» Up; L) » sin 2 29 sin 2 (^Ut- iA = sin 2 29 sin 2 U ±-\ , (2) 



where 



M w2 , 48 W] Km (3) 

Am 2 Am 2 eV 2 



is the distance at which A$ becomes 2ir. In this paper, wc shall refer to Eq.(2) as the standard 
oscillation probability and to Eq.(3) as the standard oscillation length. As an aside, we note that the 
assumptions in Eq.(l) have formally cost us the Lorentz invariance of the standard oscillation formula. 
As written, it is no longer valid in all Lorentz frames. 

The historical development of the calculation of the phase factor A$ is interesting and a lit- 
tle mysterious. The first prediction of the oscillation probability was given in 1969 by Gribov and 
Pontecorvo [20] 



=WI 2 



sin 2 26> \ sin 2 29 ( Am' 
2 ~2p~ t 
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The important point to note here is a factor two difference with respect to the relative phase which 
appears in Eq.(l). Only some years later, to be more precise in 1976, was the standard oscillation 
probability obtained by Fritsch and Minkowsky [21]. 

For a long time thereafter, the oscillation probability (2) stood as the fundamental starting point 
in neutrino physics and no comment was ever made on the additional factor two which appeared 
in the Gribov and Pontecorvo work (for further details see ref. [10]). In 1995, Lipkin re-discussed 
the derivation of the oscillation probability and pointed out that, by assuming the equal momentum 
scenario, an extra factor of two appears in the oscillation phase [4] . The hypothesis of equal energies 
was then suggested to re-obtain and justify the standard result [5]. However, simply by following 
the reasoning of Lipkin [4], the authors of ref. [8] showed that, contrary to Lipkin's assertion, the 
only case in which the standard oscillation phase can be reproduced is in the equal velocity scenario. 
This condition also distinguishes itself from the others (equal energy or equal momentum) by being 
Lorentz invariant. This docs not mean that this scenario, which yields the maximum simplicity and 
recovers the standard oscillation probability, coincides with the real situation in neutrino production. 
As observed in refs. [7,8] common velocities imply E 1 /E 2 = Pi/p 2 = m 1 /m 2 and since this may be very 
far from unity it would contradict the estimates of the neutrino energies for the known production 
mechanisms. A discussion of the kincmatical constraints derived from energy-momentum conservation 
in neutrino production will be given in Section IV. 

In order to understand how extra factors appear in the oscillation formula, let us analyze the 
difference of phase responsible for the oscillation phenomenon. In the plane wave formalism the ap- 
propriate plane wave phase is associated with each mass eigenstate. Since the four-momentum of 
different masses cannot coincide, the phase of each mass eigenstate will change with time and dis- 
tance. Thus an initially pure flavor-eigcnstate will be modified with time. The mass-eigenstate phase 
difference is 

A$ = A(ET-pL) . (4) 

In the standard treatment, one evaluates this by setting AT — AL = [22,23]. Nevertheless, if the 
two mass eigenstates have different speeds, by assuming instantaneous creation, we should experience, 
at the common time T, the interference between wave function points which have travelled different 
distances, i.e. 

L 1 =v 1 T and L 2 = v 2 T [ => AL = T Av and T = L/v ] . 

The interference between wave function points which have travelled different distances is the source 
of an extra multiplicative factor, 



_ A$(AL ^ 0) 
E ~ A$(A£ = 0) 



_ A$(Au ^ o) 
= A$(Au = 0) 



(5) 



in the oscillation phase (length). In order to quantify this effect, we explicitly calculate the difference 
of phase given in Eq.(4). By simple algebraic manipulations, we obtain 

A$ = ( -AE- Ap- ^Av ) L . (6) 

\ v v J 

As particular cases, we can immediately compute the Lorentz invariant factor A<f> in the common 
velocity and common energy scenario, 

Am 2 — Am 2 — 

A$ [ Av = ] = L w — L [standard result] , 

1 1 2p 2E 1 1 

Am 2 - Am 2 - 

A$ [AE = 0] = L « L [extra factor two] . 

p E 1 

In order to get a more general expression for the extra factor, we rewrite e in terms of p, E, Ap and 
AE as follows 

e(P,E;A P ,AE)=l + ^^. (7) 
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Considering Ap <C E and AE -C E, we obtain 



AEAp /AE 
E 2 ' \~E~ 



e(p,E;A P «E,AE«E) = i + |^Z^ |l + Q 

„ t , p(^Ap-pA^) fg) 
E (pAp - EAE) K ' 

1 when Ap = and p 12 <^ E 12 , 

2 when Ap = and p 1>2 m E 1>2 ■ 

There are difference in the value of e between the scenarios of common momentum and common 
energy but they are only significant in the non-relativistic limit. Such a situation could in principle be 
tested, for example, in the neutral kaon system. In neutrino physics, where non-relativistic neutrinos 
are unobservable, when Aw ^ 0, we practically always find an extra factor two in the oscillation length. 
In conclusion, to recover the standard formula for neutrinos in this formalism (where the energy of 
production is approximately known and is orders of magnitude greater than the postulated masses), 
in addition to the exact common velocity scenario we would also need to impose almost equal masses 
to guarantee E 1 /E 2 ~ 1 [8]. 

Now it is important to observe that, in all the above, one has implicitly assumed that the flavor 
cigcnstate is always given by the mixing matrix (chosen real by convention) with which we started 
this Section. That is, the flavor eigenstate has been assumed identical at all points and/or times of 
creation. However, we can of course multiply a flavor (mass) eigenstate by a phase factor without 
modifying its flavor (mass). Perhaps less obvious, this phase may even be space-time dependent. A 
significant example of this occurs in the next Section and is generalized in Section V. An alternative 
wave packet example, specifically devised to approximate at a given time the above non standard 
oscillating phase, is presented in Section VI. 



III. WAVE PACKET FORMALISM WITH INSTANTANEOUS CREATION 



In the preceding section, we introduced the fundamental arguments leading to extra factors in the 
oscillation probability. In this section, we are going to show why these factors do not appear in the 
usual wave packet formalism. We begin by trying to understand qualitatively the problem in a very 
simple case, that is Ap = [24,25], and deduce from it several important conclusions. 

So far, we have only considered a single plane wave. Rigorously, such an energy-momentum eigen- 
state cannot represent a physical state - it is not normalizable. It would also pose us with the problem 
of defining L and T in the oscillation phase. To avoid these problems, we must employ a normalized 
superposition of plane waves 

exp {— i [E(p, m„)t — px\\ 
and describe the time evolution of flavor states by the wave packet 

^(x, t) = ipi_(x, t) cos v 1 + ip 2 (x, t) sin# v 2 

= [ip 1 (x,t) cos 2 6 + %p 2 (x, t) sin 2 0] v a + [ip 2 (x,t) 
= ip a (x,t;6) v a + ipp(x,t;9) u , 



■01 (x, t) ] cos# sin# 



(9) 



^n(x,t) = -h= i tp n (p) e-xp{-i[E(p,m n )t-px]} dp . 



where 



As a model assumption, we suppose that the momentum distributions <p n {p) are given by Gaussian 
functions peaked around the mass eigenstate momenta p„, i.e. 



<Pn(p) = (jfaj CX P ~°"(P~P™) 2 / 4 



(10) 
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To instantaneously create at t = 0, in a localized region centered around the spatial coordinate x = 0, 
a flavor state i/ al we have to impose the following constraint 



&(a:,0) = Mx,0) [=>Aa=Ap = 0]. 
Consequently, from Eq.(9) we get 

i 

3 



*(cc,0) 



2 

7ra 2 



cxp 



.r- 



exp [ip x] u a = ip(x, 0) u a 



(11) 



(12) 



where p 

— Pi — P2- The probability P(v a — > v ; t) of observing a flavor state u p at the instant i is 
equal to the integrated squared modulus of the u p coefficient in Eq.(9), 



P(u a -> u ; t) = / \ipf3(x,t 



' dx 



26 



1 — Re / ipi{x, t) ijj* (x, t) dx 



/2 



(13) 



Actually, this result with a unique time assumes that the detector is not localized in a region smaller 
or comparable to the size of the wave packet. Otherwise for t we would have to use the average time 
of measurement. In order to calculate the oscillation probability, let us change the x-integration into 
p- integration 



t) tp* (x, t) dx 



and use the following approximation 
E(p, ni 1 ) — E(p, m 2 ) 



ip 2 {p) exp{-i[E(p,m 1 ) - E(p,m 2 )] t} dp 



Et 

p 2 -_pI 

2E 2 
(P - Po) Po 



E 1 - 
AE 
AE 



P -Po 

El 



E 2 



E 2 

w AE + (p - p a )Av . 

This approximation is justified if we assume Sp <^ E and AE <C E. The oscillation term is then given 

by 



ip*(x,t) ip 3 (x,t) dx 



exp \—i AEt] exp 



exp [— iAEt] j exp [—a 2 a 2 /2] exp [— z cr Av t] da 
Avt s " 



iV2 



By observing that AE = Am 2 /2E and using the approximation T = L/v ~ L (where T stands for 
the observation time), Eq.(13) reduces to 



P (y a — > u ; L) « sin 2 29 < 1 — cxp 



Aw T 
aV2 



cos 



A™ 2 
2E 



/2 



(14) 



Thus, when minimal slippage occurs (Av T <C a) the standard oscillation probability (2) is a good 
approximation and no extra factor appears in the oscillation term . This contradicts the result given 
in the previous Section, where, in the equal momentum scenario by using the plane wave derivation, 
an extra factor of two was obtained, see Eq.(8). To explain this apparent paradox, we observe that, 
at time T and at a fixed position x f in the overlapping region, we experience the interference between 
space points whose separation at creation is given (see Fig. 1) by 



Ax,„ = —AvT 



6 



Stefano De Leo et al: WAVE PACKETS AND QUANTUM OSCILLATIONS 



This implies that an additional initial phase, 

A$ to = -p Ax in = p Av T , (15) 

is automatically included in the wave packet formalism. Consequently, the final result contains both 
the phase difference calculated in the previous Section, i.e. 

A$ = T AE — p AL , 

and the additional term given in Eq.(15). Thus, the standard result, 

A$ st = A$ in + A$ = T AE , (16) 

is obtained. Hence, the difference in this scenario with that of the previous Section is that here the 
flavor eigenstate is not unique at all points of creation. Each point is associated with an appropriate 
x-dependent phase. 

Before proceeding further, we must ask whether the above equal momentum scenario is physically 
possible, and if not, how it is to be modified while maintaining the creation of only a flavor eigenstate. 
These modification could well change the oscillation phase. To respond to these questions we must 
first review the kinematics of particle creation. 



IV. KINEMATIC CONSTRAINTS IN PRODUCTION 

We start by observing that any production process of a particle (be it an oscillating particle or not) 
can be considered, for kinematic purposes, as an effective two-body decay such as 

M 

m < — • — ► m lj2 

where we recall that the subscript in to 1j2 refers to the oscillating mass eigenvalues. If the production 
process is a decay into more than two particles, then m represents the effective mass of all the accom- 
panying particles and is of course greater than or equal to the sum of their masses. For production 
processes other than decays M is simply the center of mass energy and not the mass of a resonance. 
In this "rest" frame, energy and momentum conservation imply 

M = y/pl + m 2 + ^Jpl + m 2 n = y/pl + ml + E n , n = l,2. 

By simple algebraic manipulations, we find [26,27] 



AE=^, Ap^^AE and A, = l^MRjlMH 
2M p p 

The next step is to observ that by assumption 



AE 
~2E 



AE . (17) 



Am ^ and thus 2(M - E) = yjp\ + m 2 + yjp 2 2 + m'l ^ . 

Consequently, in the rest frame of the decaying particle of mass M (or, in general, in the center of 
mass frame) we have 

AE ^ , Ap ^ and Av ^ . (18) 

This implies that there does not exist any frame in which Av = since this is a Lorentz invariant 
condition. We can also show that there does not exist any frame in which Ap = 0. In fact, by 
performing a Lorentz transformation with velocity (3 from the rest frame of the decaying particle M 
(or center of mass), we find 

AE' = 7 (AE - (3Ap) , 
Ap = 7 (Ap - f3 AE) . 

To satisfy Ap' = 0, we have to impose the following unphysical condition on (3 

= \Ap/AE\ = \{E-M)/p\ = {y/ft+ml + y/fi + ml) /2p > 1 . 
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This shows that (AE, Ap) is space-like. Therefore, there will, on the contrary, always be frames in 
which AE' = 0. 

It is now important to realize that the condition Ap ^ automatically implies for x ^ that the 
mass eigenstates defined in the previous Section are no-longer equal at t = 0, 

^(z,o) ¥= Mx,o) . 

This inevitably leads, within the context of instantaneous creation, to a non zero probability to find 
(see Fig. 2) a flavor state v at time t = [9, 19]. Indeed, for Ap ^ and with instantaneous creation, 
we obtain the following oscillation probability 

P (v a -> u ; t) ps sin 2 20 jl - exp 

Thus, there does not exist any time for which the state is a pure flavor eigenstate. In the next Section, 
we shall describe how by generalizing to non instantaneous creation we can eliminate the initial 
difference of phase between the mass eigenstates and hence achieve pure flavor creation event-wise . 



V. NON INSTANTANEOUS CREATION 

We have identified the initial difference of phases in the mass eigenstates 

-Apx c , 

where x c is a generic space point in the creation wave packet, as the cause for having at the time of 
creation a state which is not a pure flavor eigenstate. This undesired effect can be removed either by 
the unphysical assumption of equal momenta or by introducing for each space point x a a corresponding 
creation time t c which satisfies the following relation 

AEt c — Apx c = . (20) 

This condition guarantees 

dP(v a -> Vp\ x c , t c ) = 

and consequently allows for pure flavor creation event-wise. 

Somewhat surprisingly this substantial modification does not invalidate the standard oscillation 
formula. In fact, by following the plane wave phase calculation of Section II, we find at the interference 
space-time point (x f ,T) the following mass eigenstate phases 

=E 1 (T-t hc )-p 1 v 1 (T-t hc ) , 
$ 2 = E 2 (T~t^)-p 2 v 2 (T-t 2 , c ) . 

The last two terms in the difference of phase 

A$ = T [AE - A(pv)} - A (Et c ) + A (pvt c ) (21) 

represent the generalization of Eq. (6) in the case of non instantaneous creation. However, as explained 
in the previous Section, in the wave packet formalism additional initial phases arc automatically 
included in the expression of the oscillation phase. Thus, for non instantaneous creation, we still have 
to take into account the contributions of the initial phases 

= £A,c - Pi [x f —v 1 (T— t lje )] , 
$2,in = E 2 t 2 , a - p 2 [x f —v 2 (T— t 2)C )] . 

The final result contains both the difference of phase (21) and the additional term 



A$,„ = -Apx f +TA(pv) + A (Et c ) - A (pvt c ) 



(22) 
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Finally, after integration (x f — > L) the standard result 

A$ st = TAE-LAp (23) 

is once more obtained. 

The above procedure has eliminated flavor contamination at creation. Nevertheless, since creation 
no longer occurs at a unique time and the partially formed wave packets naturally evolve, there will 
still not be a pure flavor eigenstate at any fixed time, with the trivial exception of the very first 
instant in the creation process. We also observe that any search for the particle during this time 
(creation) will not necessarily yield a positive result since the wave function is not fully normalized. 
The measurement will still produce a collapse of the wave function in the appropriate percentage of 
cases to zero. 



VI. GEDANKEN WAVE PACKETS 



We have seen that the implementation of pure flavor creation, while non trivial, does not modify the 
standard oscillation formula. However, the wave packet assumed was by no means the most general. 
Now we shall study the possible consequences of a two-peak wave packet. This is the simplest which 
allows for the insertion of additional constant initial phase factors. 

To simplify the following calculation, we return to the unphysical (in any frame) Ap = scenario 
with instantaneous creation at t = 0. We consider a wave packet obtained by a sum of generalized 
Gaussians with peaks at x = ±x (see Fig. 3), 



tp(x,0) = N < exp [ip (x + x )] exp 



(a; + x ) 2 



- exp [ip (x — x )] exp 



x y 



(24) 



where N is the normalization constant. We also assume that the peaks are well separated. In this case 

1 



N 



2-ko? 



Note that each gaussian has its own extra constant phase factor. 

We now allow the wave packets to evolve and consider the situation (measurement) in which the 
first peak of one mass eigenstate overlaps with the second peak of the other (see Fig. 4) . In this region 
of overlap, oscillation occurs but it is to be noted that there are non-overlapping parts of the wave 
packet which correspond to pure mass eigenvalues (see Fig. 5). The contribution of these parts yields a 
non oscillating term. Indeed this is exactly what happens when the mass wave packets have completely 
separated (decohcrence) . The oscillating term is modified by the presence of the difference between 
the phase factors exp[ip x ] and exp[—ip x ]. Thus, the oscillation formula at or around these times 
reads 



i sin 2 29 \ 1 - \ exp 



Avt- 2x 
aV2 



Am 2 
2E 



/2, (25) 



where the sign in the constant term =p 2p x depends on v 1 ^v 2 . In this picture, 

2p x fn ± pa Av t « =F -r-pr t . 

Thus, we obtain the extra factor two of Section II for the oscillating phase. It is true that since 
T-2p x is a constant term the oscillation period is still be the standard one (see Fig. 5). However, 
with multiple peaks, each with independent initial constant phases, the effective extra phase term 
could become L or T dependent. We see no good physical reasons for excluding such a contribution 
a priori. 

In the case of minimal slippage when the interference occurs within each separate peak, the constant 
phases play no role. In these cases, we can also ignore (by definition) the non-overlapping parts of the 
wave function. Thus, the standard result is obtained as a good approximation. This is now consistent 
with the fact that for no slippage (Ad = 0) the standard formula is exact. 
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It is also useful to observe that, in the incoherent limit, one has a clean method for determining 
experimentally the mixing angle, whatever the wave packet form is, 



Nevertheless, the time for the onset of complete incoherence depends upon the wave packet dimensions 
and upon the different mass velocities. 



The aim of this work was to shed some new light on the quantum mechanics of particle oscillations. 
The primary objective was to search for the conditions under which the standard oscillation formula 
is valid. In the process we have understood the origin of the extra factors in the plane wave oscillation 
calculations. It is the implicit assumption that at creation (whether instantaneous or not) the flavor 
cigenstate is unique even up to the phase at all points and times of creation. To the best of our 
knowledge this has not been pointed out previously. The often quoted plane wave derivations of the 
standard formula have generally been based upon invalid approximations or formalisms chosen so as 
to compensate for the neglect of the initial phase contributions. 

Other authors have previously pointed out that the assumptions of equal velocity, momentum 
and energy are "unphysical" in the sense that they are not compatible with the known production 
processes in the laboratory frame. We have pointed out here that the first two are rigourously non 
physical in the sense that there is no Lorentz frame in which they occur. Only the equal energy case 
is theoretically possible. Consequently, since the assumption of instantaneous creation together with 
flavor eigenstate production imposes equal momentum, it is non physical. We can correct for this by 
assuming an event-wise production mechanism. Event-wise production is perfectly natural. It is even 
predicted for the equal momentum case (had it existed) when seen by another observer in a Lorentz 
frame in which the momenta are not equal. Somewhat surprisingly the standard oscillation phase, as 
calculated in the wave packet formalism, is not affected by this modification, at least not for the cases 
with minimal slippage. Thus, we have concluded that within the wave packet formalism the standard 
result is not only exact in the case of equal velocities (no slippage) but also a good approximation in 
all cases in which minimal slippage occurs between the mass wave packets. 

Now the standard oscillation formula contains no dependence upon the form of the wave packets 
involved. This clearly cannot be valid in general. Indeed, as a simple and well known counter example 
we have recalled in the previous Section the incoherence limit, who's onset is dependent upon the 
size of the wave packets and the differences in the mass eigenstate velocities. One must expect the 
dimensions and explicit form of the wave packets, including all relative phases, to play a role in the 
oscillation formula. We have exhibited a simple two-gaussian model to demonstrate not only this 
but also how a result that simulates the extra-factors calculations (for the oscillating part) can be 
obtained. 

Allowing for completely arbitrary modulations of the plane waves, we must introduce additional 
parameters into our generalized formula 



where R(t) and <p(t) depend upon the details of the wave packet envelope. This formula contains two 
extra parameters when compared to the standard formula, four in all. It may well prove necessary to 
employ our generalized expression in order to reconcile diverse experimental results. The alternative 
could lead to inconsistencies in the determination of mass differences and/or mixing angles. Use of this 
formula may also avoid the need to introduce one or more sterile neutrinos. This equation, as it stands, 
is probably too vague for practical phcnomenological fits. However, it is always possible, without 
returning to the standard formula, to add simplifying assumptions such as the time independence of 
the extra phase term. 

Apart from the question posed in this paper of what is the most practical oscillation formula 
to be used in phcnomenological fits, there is an aspect of this work which is of great interest, at 
least to the authors. The details of the creation and annihilation of wave packets is to a large degree 
unknown territory. Oscillation phenomena may indeed be useful as a source of information upon this 



v ; T>a/|Av|) = sin 2 29/2 . 



VII. CONCLUSIONS 




(26) 
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subject. It should also be possible to investigate this aspect for photon wave packets with the help of 
very precise measurements in interference phenomena. In the case of photons the effect of slippage is 
substituted by the occurrence of different path lengths. For example, in optics it is well known that 
interference effects cease if the difference of path lengths exceeds the wave packet dimension, and this 
should permit the determination of these dimensions. This technique may, of course, also be extended 
to any elementary particle that lives long enough. This particular subject matter recalls transitory 
phenomena in various sectors of classical physics. Its study, both theoretical and experimental, at the 
quantum mechanical level is surely a great challenge. 
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Fig. 1. The square modulus of the mass eigenstate coefficients, \ipi(x,t)\ 2 (upper-half) and \ip 2 (x,t)\ 2 (lower- 
half), is plotted as a fucntion of x for two times: t = (left-side) and t = T (right-side). As a model assumption, 
the wave packets are supposed to be Gaussian functions (with the same width a) peaked around x — 0. At 
observation time T, the mass eigenstate wave packets are centered around different space-points, Li — ViT 
and 1/2 = v 2 T (spreading effects are neglected). Since Li and L 2 are assumed very large compared to a, the 
plots of the mass eigenstate wave packets at creation and observation are separated by dots in the x-axis. At 
observation time T, Xj is a fixed point in the overllaping region. At this point, the mass eigenstate phases are 
given by $ 1 (x f ,T) = E 1 T—p 1 v 1 T and $ 2 (x/,T) = E 2 T — p 2 v 2 T. The interfering wave packet points at (x f ,T) 
correspond to different initial wave packet points. Thus, the initial phases &i,i n (%j — ViT, 0) = —piX/ +P1V1T 
and <&2,in{xf — v 2 T, 0) = —p 2 Xf+p 2 v 2 T are automatically included in the wave packet formalism. Consequently, 
the standard result TAE — x f Ap is obtained. 
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